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Abstract 

A canonical ensemble model is used to describe a caloric curve of nuclear 
liquid-gas phase transition. Allowing a discontinuity in the freeze out density 
from one spinodal density to another for a given initial temperature, the 
nuclear liquid-gas phase transition can be described as first order. Averaging 
over various freeze out densities of all the possible initial temperatures for a 
given total reaction energy, the first order characteristics of liquid-gas phase 
transition is smeared out to a smooth transition. Two experiments, one at low 
beam energy and one at high beam energy show different caloric behaviors 
and are discussed. 



1 



Typeset using REVTgX 



Recent experiments of heavy ion collisions at low and high beam energy show caloric 



also used in extracting the critical exponents of nuclear system [Q. 

The investigations of nuclear phase transition through heavy ion collision are based on 
two basic assumptions which are not confirmed: One is that we can apply equilibrium ther- 
modynamics for such a small system of few hundreds constituents at the most. Another is 
that a thermalized uniform system is formed in heavy ion collision before the multigfrag- 
mentation takes place. Based on these assumptions, a canonical ensemble model describing 
nuclear multifragmentation phenomena have been developed ||. Another assumption of 
this model is that the fragment distribution is directly related with the state of thermalized 
uniform system before its break up. In this paper we will use this model to describe the 
caloric curve of nuclear liquid-gas phase transition. This model may also determine a beam 
energy required for the study of critical behavior. 

We first summarize some results of the canonical ensemble model of Ref. || which we 
will use in this paper. Further details and derivations can be found in that reference and 
references therein. In Ref. ||, when each cluster of size k is weighted by Xk = iXf, the 
canonical ensemble averaged number of clusters of size k is given by 



for a system of A nucleons. The result is independent of X±, i.e., the weighting factor 
assigned to each nucleon. The mean cluster distribution thus depends only on x which is 
the same for each cluster. As x — > 0, the system forms one large cluster. As x —>■ oo, only 
nucleons are present. 

The Xk can be obtained || by weighting each cluster with its free energy as Xk = 



curves for nuclear liquid-gas phase transition that are different Jl|J3]. Even if there is some 
uncertainty concerning the extraction of temperatures in heavy ion collision ||, low energy 
collision shows some indication of a sharp first order phase transition |IJ] while higher energy 
collision shows a much smoother phase transition |2j . The data of higher energy collision is 
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Here Mb is the nucleon mass, g = 4 is the degeneracy factor. The freeze out density p is 
the density when there are no further break up or clusterization occuring. An equilibrium 
temperature T determines the energy shared between the liquid phase part (internal energy 
of clusters) and the gas phase part (thermal motion of clusters). If we approximate the 
energy for a cluster of size k at temperature T as 
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then the common weight factor of each cluster becomes x = Ax p T 3 / 2 . Here is the average 
binding energy at zero T, eo is the level spacing parameter and the T is a cut-off temperature 
for internal excitation of the cluster. The a# ~ 16 MeV for a nuclear matter and a# ~ eo ~ 8 
MeV for finite nuclei. 

Substracting the thermal energy |T of the center of mass, the canonical ensemble aver- 
aged total energy per nucleon is given by 
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The excitation energy per particle E*(T) = E{T) + as is 
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The low T limit of Eq.([5D, E(T —> 0) ~ T 2 /eo, gives the Fermi liquid limit and the high T 
limit of Eq.flD, E{T -> oo) w §T, gives the ideal gas limit These show that the caloric 
curve is independent of the freeze out density for low and high temperatures. 

Eq.flD) shows that the phase transition diagram (E-T caloric curve) of a nuclear sys- 
tem is determined once the freeze out density p is known. Eq.(|^) is a smooth function of 
temperature T for fixed p with a constant x p . If the nuclear density p changes smoothly as 
the T changes, then E(T) will have an extra T dependence than the explicit dependence 
of Eq.(^[) but it still changes continuously. However, if the density p changes suddenly at 



some temperature T, then the excitation energy Eq.(]5|) has a discontinuous jump at the 
temperature through the discontinuity of x p . The gap of this jump is the latent heat of 
the liquid-gas phase transition. In a spinodal instability region, a system can not keep its 
uniformity and thus breaks up into a mixture of gas and liquid. The densities at the liquid 
and gas spinodal points of an isothermal curve for a given T gives a jump in the freeze out 
density from liquid to gas phase. 

For a thermally equilibrated uniform system, the density is determined once we know 
the temperature and pressure or the temperature and energy as shown in Fig.[l]. If we add 
energy to the system keeping the temperature fixed, then the density and the pressure will 
be changed to new values according to the isothermal curve of Fig.[l| On the other hand, if 
we add energy keeping the pressure (density) fixed, then the temperature and the density 
(pressure) will be changed. To maintain non-zero pressure for a finite system, we need some 
mechanism which can controll the volume such as infinite potential wall. However, since such 
a mechanism does not exists, the system should move toward the P = point accompanied 
by a corresponding temperature and density change. Thus if we add energy adiabatically 
to a finite nucleus, then the system changes along an isobaric path with zero pressure. Such 
a system with conserved total energy E will finally settled down with some temperature T 
and zero pressure P. This process is similar to the equilibration of (adiabatically) heated 
water at atmospheric pressure. 

As shown in Fig.|| the pressure is always positive for T > T z where T z is the temperature 
at which the minimum pressure is zero. In the lower figure of Fig.|], the P = isobaric 
curve is shown by a dashed line. Nuclear matter is always in spinodal stability region for 
T > T c where T c is the critical temperature. The spinodal curve is shown by dotted line in 
this figure. Also shown in Fig.|TJ are the critical point c, the point n of maximum energy E n 
of spinodal curve, the point m of maximum energy E m of zero pressure isobaric curve, and 
the point z with energy E z which is a spinodal point with P = 0. 

When nuclear matter is heated adiabatically, the caloric curve should follow the zero 
pressure isobaric curve shown by dash-dotted line in Fig|| until it reaches the point of 
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T = T z where the zero pressure state gets into the spinodal instability region, i.e., while 
E < E z . For Fig.|I|, where a B = 15.77 MeV, T z = 15.16 MeV and E* = 17.7 MeV. If 
the energy E becomes higher than E z , the system can not be in a liquid phase as a whole 
with P = 0. The only possibilities are then a uniform gas phase of constituent particles 
or a mixture of clusters of liquid phase and constituent particles. If a thermalized system 
belongs to the spinodal instability region, then the system would fragment into clusters. 
For dynamical equilibrium (stability) of a cluster, each cluster should have zero internal 
pressure. For thermal equilibrium, the T associated with thermal motion of any cluster and 
its internal excitation should be the same. For a thermally and dynamically equilibrated 
system, the mean fragment distribution and the ensemble averaged total energy are given 
by Eqs.dH) and (|5|) respectivley. For E > E n , the system should expand to infinite volume 
approaching an ideal gas phase of nucleons since it stays outside the spinodal instability 
region with positive pressure. 

A system produced in heavy-ion collision can break up into pieces before it reaches 
thermal equilibrium. During such a collision, the nucleons and the clusters broken off take 
some energy as separation energy and thermal energy and the remnant equilibrates thermally 
at lower T and E. For a heavy ion collision with beam energy per nucleon of Eb eam , the 
thermalized remnant would have E = E* — a# and T within the limits < E* < Eb eam and 
< T < T max , where T max is the maximum T which can be reached by a uniform system 
with energy E, i.e., the minimum energy of the T max isothermal curve is E. 

Suppose we have a thermalized uniform system with energy E and temperature Tj sat- 
isfying EOS of Fig.|l|. For a given value of E and T i} there are two possible values for p^ and 
Pi unless E > ^Ti which has only one possible value at pi > p . For a given temperature T i; 
which is lower than the critical temperature T c , there are two spinodal points; one at the 
liquid side with pi and E\ and another one at the gas side with p g and E g . We now discuss, 
assuming equilibrium thermodynamics, different possible behaviors of the remnant system 
in various energy regions for a given initial temperature Tj. 

If E is low enough to be in the liquid side of the spinodal stability region of the Tj 
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isothermal line, i.e., E < Ei, then the system is either compressed or expanded (due to 
its nonzero pressure). However, the system maintains uniformity since it is in the spinodal 
stability region and the temperature should change to keep E conserved. For E < Ei, the 
system will eventually equilibrate to the zero pressure liquid phase at Tf with energy E. If 
E = Ei with pi, the system equilibrates at Tf = 7]. Here TJ is the temperature at which the 
energy of the P = point is E\ and it is obvious that Tj < 7] < T max (see Figjj]). In this 
case no break up has occured and the system equilibrates at temperature 7]. 

If now E is such that E\ < E < E g , then the system is either in the spinodal instability 
region of Tj or is in a highly compressed dense liquid phase. If the system stays uniform long 
enough, then it will be compressed or expanded depending on the pressure until it reaches 
Tf with P = and energy E. For this situation, Tf > 7). However, for < p , the system 
would not have a uniform density due to the spinodal instability. Instead, the system will 
break up, thereby becoming non-uniform and composed of two parts with the total energy 
E conserved. One part will be a liquid phase with E\ and p% and another part is a gas phase 
with Eg and p g at temperature Tj. We may suppose that the system breaks up (fragments) 
at p < pi and the constituents clusterize at p > p g which becomes the freeze out densities. 
The energy may flow from the dense part to the dilute part, appearing as liberation energy 
and the thermal energy of the fragments during this dynamical fragmentation process. If we 
assume that the dynamical rearrangement process is faster than the thermalization process 
in the spinodal instability region, then the temperature stays at T while it fragments. 
However, the clusters of liquid phase will eventually equilibrate internally to Ei, T\ and 
zero pressure without any further clusterization or fragmentation. The gaseous system 
which is composed of nucleons and clusters (each cluster behaving as a particle) would then 
equilibrate thermally with the internal state of clusters at 7}, and then it expands to an 
infinite volume due to the positive pressure. The 7} becomes the freeze out temperature of 
the gas. Here Tf = T\ and E = aEi + (1 — a)E g but pi ^ api + (1 — ot)p g . 

If the energy is even higher such that E g < E < E n , then the system is either in the 
uniform gas phase with density pi lower than p g or in highly compressed very dense liquid 
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phase with positive pressure. If pi < p g , the system expands to an infinite volume to achieve 
zero pressure without any clusterization. If the system were initially in a very dense liquid 
phase pi > po, then it can reach a zero pressure point with temperature Tf within the 
liquid phase of the whole system (for E < E z ). Another possibility is an expansion from 
the dense phase that over shoots into the spinodal instability region or to the gas phase 
region of the temperature %. If it ends up in the gas phase region, i.e., p < p g , then it will 
further expand to infinite volume since P > 0. If the original expansion reaches the spinodal 
instability region, then the system would break up into a mixture of liquid phase and gas 
phase. However since the energy E is higher than E g , the liquid phase should share more 
energy than E\ to thermally equilibrate the internal temperature of each cluster of liquid 
phase and the temperature of thermal motion of the gaseous fragments. Thus the freeze 
out temperature Tf now become higher than and Tf increases as E increases further for 
E > Eg. Here E > aE\ + (1 — a)E g and the exceeded energy of this inequality makes 
T f > 2J. 

If the energy is higher than E n , then the system would not have any region of spinodal 
instability. Thus the only possible phase is the gas phase (due to P > 0) made of nucleons 
without any cluster of liquid phase. 

These processes are represented by each broken line for each temperature Tj in Fig.|2|. If 
Pi > pi with E < Ei, then the whole system becomes a uniform liquid at Tf having zero 
pressure and energy E. The low energy - low temperature part of Fig.|2] represents the zero 
pressure point of a liquid phase which is very close to E* = T 2 /eo, the low T limit of Eq.(^). 
For the nuclear matter equation of state shown in Fig.|l|, E* = 2.1 MeV (point s) and the 
T\ pa 5 MeV for Tj = 0. The liquid side spinodal density p\ becomes the lower limit of the 
spinodal stability region of the equilibrated Fermi liquid system and T} ~ 5 MeV becomes the 
break up temperature H at which a nuclear system fragments into two clusters in average. 
The thermal properties of a nuclear system below E ~ 2 MeV and T pa 5 MeV are hard 
to be extracted from data |2J. If pi < p g with E > E g , then the whole system becomes a 
uniform gas at the freeze out temperature Tf and density p g . The very high temperature - 
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high energy part (which is not shown in Fig.Q) represents a gas system without any clusters 
and the high T limit of Eq. (|^) applies to this region. 

If the density is p g < pi < p\ then the system should be a mixture of liquid of p > p\ 
and gas of p < p g at temperature Tf with zero pressure. The spinodal densities pi and p g of 
temperature T become the break up (fragmentation) density of clusters and the freeze out 
(clusterization) density of gas respectively at Tf. For E < E g , there is not enough energy 
for the system to become a uniform gas, and for E > Ei, there is too much energy for the 
system to maintain a uniform liquid phase as a whole. The flat part of the caloric curve of 
Fig.||] represents the spinodal instability region with E\ < E < E g between pi and p g where 
some parts form clusters of liquid phase of Tf = Ti with E\. In turn, this system of clusters 
forms a gas system at 7], with freeze out density p g . The higher temperature - higher energy 
part of Fig.|] represents the liquid-gas mixed phase of temperature Tf > T\ with E > E g . 
Here the liquid phase has zero pressure with energy higher than E\ and the gas phase has 
energy E g and freeze out density of p g . For the nuclear matter results of Fig.[TJ, E* = 16.15 
MeV for zero temperature. This behavior of E-T caloric curve can be represented by the 
results Eq. (|5|) with two different freeze out densities, p\ for liquid phase and p g for gas phase. 
As can be seen from Fig.^], the data of Ref. [|TJ can be fit with Tj = ~ 1.5 MeV and with 
a B = e = 10 MeV. 

For a given energy E, the temperature T of system may have any value consistent with 
< T < T max . Various T cause various Tf and thus we need to superpose these various 
temperatures with some weight. For smaller beam energy the temperature range is smaller 
and the average is not much different than taking a single temperature. A superposition 
of T's with < T < 2 MeV and with a B = e = 10 MeV fits the data of Ref. of 600 
MeV/ A Au + Au collisions well (the thin solid line which is similar to the thick solid line 
for T = 1 MeV using the spinodal densities 0.677p and 0.015p as the freeze out densities 
of liquid and gas phases respectively). For a collision with larger energy, averaging over a 
larger range of temperature gives a smooth curve of E-T. The superposition < T < 8 
MeV with a B = eo = 6 MeV fits the data of Ref. fH of 1 GeV Au + C collisions (the smooth 
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thin solid line). This range of T\ is compatible with T{ = 8.7 MeV at critical multiplicity 
given Ref. 0. 

The above analysis indicates that the liquid to gas phase transition (break up of nucleus) 
starts to occure at the temperature and energy at which the pressure is zero and the energy 
is the same as the liquid side spinodal point of zero temperature. It also indicates that a low 
energy collision would show a clearer first order phase transition than a high energy collision 
once the excitation energy is high enough to reach the zero temperature spinodal energy. A 
high energy collision with T max > T c is required in study of critical behavior. The different 
set of values of a# and eo used in the fits for the low and the high beam energies show 
that the multifragmentation process with lower beam energy is better applicable than the 
higher energy reaction in extracting information for infinite nuclear matter. The lower energy 
reaction gives larger size fragments and thus closer to infinite system. We need further study 
to check if the dynamical rearrangement process is faster than the thermalization process in 
the spinodal instability region and the other way around for spinodal stability region. 
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FIG. 1. Equation of state of nuclear matter || for temperature T of, from bottom to top, 
0, 3, 5, 7, 10, 15.16 (T z at which the minimum pressure is zero), and 20.95 (which is the critical 
temperature T c ) in MeV. Also shown are the spinodal points (dotted line), the points with zero 
pressure (dashed line), and the points of minimum energy (filled circles) for each temperature T. 
See the text for variously marked points of cross. 
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FIG. 2. Phase transition caloric curve. The open circles are the data of Ref. |IJ and the crosses 
arc the data of Ref. @. The solid curves are the fit of Eq. @ with spinodal densities at T of, from 
bottom to top, 0, 1, 2 MeV, and T c with as = en = 10 MeV and To = T c for A = 50. The dashed 
curves are the fit of Eq.(||) with spinodal densities at T of, from bottom to top, 0, 4, 5, 7 MeV, 
and T c with as = eo = 6 MeV and To = T c for A = 110. The thin solid curves are the fits as 
discussed in the text. The dotted curves are, from the top to bottom, T = VlOE*, T = V6E*, 
and T = (2/3)(E* — 2.0) in MeV, as in Ref. Q. The dash-dotted line is the zero pressure curve of 
Figjl] with a B = 15.77 MeV. 
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